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Mixtures of nanoparticles and polymer-like objects are encountered in many nanotechnological 
applications and biological systems. We study the behavior of grafted polymer layers decorated by 
nanoparticles that are attracted to the polymers using lattice gas based mean field theory and accom- 
panying coarse-grained Brownian dynamics simulations. We find that the presence of nanoparticles 
can induce large morphological transitions in the layer morphology. In particular, at moderate 
nanoparticle concentrations, the nanoparticles cause a reduction in the height of the polymer layer 
above the grafting surface, which occurs via a novel first-order phase transition for sufficiently strong 
attraction between the polymers and the nanoparticles and smoothly for weak attractions. The pre- 
dictions of the theory qualitatively agree with the observed behavior of grafted natively unfolded 
protein strands upon binding of proteins [T] . The results also inform ways of designing nanopolymer 
layer morphologies. 



I. INTRODUCTION 

Polymeric materials play an important role in ma- 
terials science and engineering. Along with traditional 
polymeric materials such as rubbers and plastics, thin 
films of polymers, often grafted to surfaces to form 
"polymer brushes," have become prominent in many 
technological applications such as those reviewed by 
Senaratne in Rcf. 2J. Applications include stabilizing 
a disperse colloidal suspension [31 H] or creating de- 
vices like novel electrodes [5], organic solar cells [6], 
or organic memory [7]. Overall, the properties of such 
thin polymer films have been well characterized ex- 
perimentally and reasonably well understood theoret- 
ically (e.g. [HHH]). 

Thin polymer films decorated with nanoparticles 
are a novel class of composite materials which allow 
one to design their function and control their morphol- 
ogy in nanotechnological applications such as solar 
cells and programmable memory devices El [13 [S] . 
A class of such additives, sometimes called "sticky" 
antiplasticizers, which bind to polymers, has been 
studied experimentally (e.g. [TTHTO] ). These addi- 
tives form hydrogen bonds with the polymers to which 
they are added and thus decrease the free volume 
in the material |17j . Their addition results in mod- 
ification of the crosslinking in the polymer network 
and changes in the cooperative motion between poly- 
mers [IH] . Practically speaking, they modify the crys- 
tallinity, glass transition temperature, toughness, and 
elastic moduli of polymeric materials. Only recently 
has the addition of such additives to a polymer melt 
been studied in molecular simulations [201 121] ■ 

Polymer-like biological molecules, such as natively 
unfolded proteins, play an important role in the func- 
tioning of living cells. There are many examples, such 
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as those reviewed by Uversky in Ref. [2?. Examples 
include newly-made proteins in the endoplasmic retic- 
ulum [221 , a protein implicated in Alzheimer's Disease 
called NACP [24|, and nucleoporins which facilitate 
transport through the nuclear pore complex (NPC) 
[5SH2H]- In the case of nucleoporins, it is common for 
these biopolymers to bind to other proteins [2^1 [30] , 
which we think of as analogous to nanoparticles. For 
instance, it has been hypothesized [3 |3I1 [31] that 
transport through the NPC involves a conformational 
change in the layer of natively unfolded proteins that 
line the NPC passageway, induced by the binding of 
certain transport factor proteins. Such changes have 
been observed in in vitro experiments [T]. However, 
their precise role in the transport of cargoes through 
the NPC is still not understood. 

Without specific attractive interactions between the 
nanoparticles and the polymers, they effectively re- 
pel each other, and the effect of this repulsion on 
nanoparticle aggregation is an important element in 
controlling the physical properties of polymeric mate- 
rials (see e.g. [TS] and references therein). Kim and 
O'Shaughnessy [T3] derived the spatial distribution of 
such nanoinclusions in a polymer brush — for instance 
whether they interpenetrate the brush or segregate 
into a layer on top of the polymers. Similar ideas were 
examined in computer simulations by Binder and col- 
leagues [14] . 

The case in which nanoparticles are attracted to 
grafted polymers has received less attention. Marko 
considered a polymer brush in the presence of a binary 
mixture of solvents |33j . which could, in principle, 
model this system. However, he focused primarily on 
the case of repulsive interactions and the case in which 
attractive interactions occur via higher order terms in 
the density expansion as compared with the case con- 
sidered here. More recently, in Ref. Halperin et 
al. considered, among other things, nanoparticles that 
bind weakly to a polymer brush. The weak binding 
assumption results in a behavior devoid of the many 
features described in our work, such as global mor- 
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phology changes. 

The main goal of the present study is to gain insight 
into the behavior of polymers grafted to a flat surface 
interacting with nanoparticles which can bind to the 
brush. We consider a situation in which the polymer 
layer is in equilibrium with a solution of nanoparticles 
which can penetrate into the layer, thereby changing 
properties such as the height of the layer above the 
grafting surface and its density. Our goal is to estab- 
lish how the properties of the polymeric layer, particu- 
larly its height above the grafting surface, depend on 
the number of nanoparticles bound inside the layer, 
the strength of their interaction with the polymers, 
and ultimately the concentration of the particles in 
the solution above the layer. We will model the sys- 
tem using a lattice gas based mean field theory and 
use this theory to describe the collapse of the poly- 
mer chains upon binding. We also will present Brow- 
nian dynamics simulations, which were performed in 
order to validate and further elucidate the behavior 
predicted by the mean field theory. 

II. MEAN FIELD THEORY 

In this section we present a general theory of the 
behavior of polymer layers mixed with nanoparticles. 
We will first obtain a free energy for the system using 
an approach based on a lattice gas model J35], and 
then we will examine in detail the features of this free 
energy and the equilibrium states predicted by it. 

A. Free Energy 

Consider a layer of Np polymeric chains, each con- 
taining N monomers, which are grafted at one end to 
a planar surface. In the absence of direct interactions 
(beyond excluded volume) between the monomers of 
the chains, such planar layers are known as "brushes" 
and have been extensively studied (e.g. [HHH])- 

Our initial goal is to determine the free energy of 
the polymer layer as a function of the essential variable 
quantities, namely the layer's height h and the num- 
ber of nanoparticles in the layer volume Mg. We 
will obtain the entropy by first considering a different 
but related system comprised of Nm = NpN discon- 
nected monomers and A^^ nanoparticles in a volume of 
elemental cells (each with volume where b is the 
monomer diameter), and then correcting it for the ef- 
fect of polymer connectivity (including grafting). One 
could also obtain a similar expression for the entropy 
using other methods (as discussed below). 

Here, we consider the case in which the nanoparti- 
cles have the same size as monomers. We have also 
considered the more general case in which nanoparti- 
cles are larger than monomers [35] • We find that the 
mean field theory predicts the same qualitative fea- 
tures, including the rich behavior of subsections |II B| 
and |II C[ We therefore focus on the simpler case for 
the sake of clarity in the present work. 



In a simple lattice gas model for a system of Ng 
nanoparticles and Nm disconnected monomers of the 
same size, the total number of configurations is 



The entropic contribution to the free energy, in units 
of kT, is therefore given by Fent = — In $7. Making 
use of Sterling's approximation and converting to in- 
tensive variables, 

+(l_0_^)ln[l_0_V^]} (2) 

where ip — Nm/M^ is the density (or, since we ex- 
press all lengths in units of 6, the volume fraction) 
of monomers in the polymer layer and (f) = Ni/Mi 
is the density (or volume fraction) of the nanoparti- 
cles. This expression includes contributions from the 
translational entropy of monomers and nanoparticles, 
as well as excluded volume in that no two particles 
may occupy the same lattice site. 

However, the monomers are not a gas, but rather 
are linked together into polymer chains, which, in ad- 
dition, are tethered to a flat surface. To correct the 
above results for polymer connectivity and grafting 
[HI EZIj the translational entropy of the monomers 
should be replaced by the expression for the config- 
urational entropy of an ideal polymer chain grafted 
onto a flat surface as in a polymer brush, namely 
—Nph?/{2N) [5]. Strictly speaking, the translational 
entropy of the monomers is —Mnil^ihuj}— 1), but since 
adding a constant to the free energy will not change 
the equilibrium states, we will, for simplicity, treat 
—Miiphnp as the translational entropy of an ideal 
gas of monomers. For conveniece, we will henceforth 
implicitly normalize free energies by the number of 
chains Np. The entropic contribution to the layer free 
energy per polymer chain (in units of kT and denoting 
M = Mi/Np) then becomes 

Fe„t = h'/{2N) + M{(j,\n^ 

+ (l-0-7A)ln[l-0-V]}. (3) 

It should be noted that our lattice gas based deriva- 
tion is not the only way to obtain the entropic part 
of the free energy of the system. A derivation based 
on Flory's approach [38l [39] , in which monomers are 
placed on adjacent lattice sites in accordance with 
polymer connectivity, would give a free energy sub- 
stantially the same as Eq. (|3|, with only inessential 
differences such as terms that do not contribute to 
its relevant intensive thermodynamic quantities and 
terms which disappear in the limit 3> 1. A van der 
Waals gas based approach for a multicomponent sys- 
tem 40 , which does not rely on a lattice, can also be 
used. The van der Waals analog of Eq. ([s]) would be 
slightly different in appearance, but would have sub- 
stantially the same virial expansion to second order in 
concentration of monomers and nanoparticles as well 
as similar behavior at high concentration. 
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Also note that when (/) = in Eq. ([s]) (no nanopar- 
ticles in the polymer layer) 

F,r,t^h'/{2N)+M{l-^)\n[l-i,]. (4) 

If the logarithm is expanded for <C 1, this becomes 
the standard mean field free energy of a plane-grafted 
polymer [HI [H], and it yields the standard brush 
height scaling behavior h - Na'"^/^ [HI H [12]. We 
choose not to expand the logarithm because the den- 
sity may become large once nanoparticles are added. 
The first term in Eq. Q describes the entropic cost 
of stretching the grafted polymer segment. The 
second term accounts for excluded volume between 
monomers, with a built-in upper limit to the density 
■0 = 1. (At this density monomers fill up the entire 
lattice, and thus the polymer layer cannot be com- 
pressed any further.) 

To include the contribution to the free energy due to 
the attraction between monomers and nanoparticles, 
we will write for the net interaction energy between 
one nanoparticle and the polymers x(c&)V'j where x 
is a phenomenological interaction constant which de- 
pends on the strength of the microscopic interaction 
between the monomers (discussed in subsection 
IV A[ below) . Thus, if there are Ng = M^cj) nanoparti- 
cles in the polymer layer, the total internal energy of 
all such interactions is, in units of kT and normalized 
by the number of chains, 



E = Mx(t>'il}- 



(5) 



Note that x < corresponds to attractive interac- 
tions. Finally, the total free energy of the polymer 
brush layer is Fb — E + Fe„t, with i^e„t given in 
Eq. (|3|, i.e., 



Fe 



h^/{2N) 
+(1-0- 



0) ln[l - 



)ln0 



(6) 



On one side of the polymer layer is the grafting 
plane. On the other side ("above" the layer) is a so- 
lution of nanoparticles. In addition to the free en- 
ergy of the polymer layer, we also need to include the 
free energy of the solution that is in equilibrium with 
the polymer layer. We will model the solution as an 
ideal solution comprised of nanoparticles mixed with 
an unstructured solvent. The free energy of the solu- 
tion, Fs, is standard, and is given by the free energy 
of mixing of a gas of nanoparticles, and hence by the 
r.h.s. of Eq. ^ with -0 = 0. 

Denoting Ms as the volume of the solution (i.e., the 
number of monomer sized cells that comprise the solu- 
tion volume) , and Ns as the number of nanoparticles 
in the solution, we compute for the solution chemical 
potential 

M5 - OFsiMs, Ns)/dNs = In [Co/(l - Co)] , (7) 

where Cq = Ng/Mg is the density (or, equivalently in 
our units, the volume fraction) of nanoparticles in the 



solution. Similarly, the osmotic pressure of nanopar- 
ticles in the solution is given by 



-dFs{Ms, Ns)/dMs = - In [1 - C'o] . (8) 



In the same way, the chemical potential and osmotic 
pressure in the polymer layer can be obtained from 
Eq. ^ as 



MB =XV' + ln[0/(l-0-V')] 



ijj — ln(l — (j) — 



(9) 



(10) 



In thermodynamic equilibrium, the chemical poten- 
tial of the nanoparticles and the osmotic pressure must 
be equal in the polymer layer and in solution, which 
are assumed to be macroscopic. Consequently, the 
chemical potential of the nanoparticles in the layer 
and osmotic pressure of the polymer layer are de- 
termined by the bulk solution concentration Co, as 
indicated in Eqs. ([t]) and Q above. Equating the 
chemical potentials and osmotic pressures of the two 
phases gives two equations in the two unknown vari- 
ables M and Nb (where Nb = N^/Np). Simultane- 
ously solving these two equations yields the polymer 
layer's volume and number of bound nanoparticles at 
the extrema of the free energy. 

One can ensure that the solution found is the global 
minimum of the free energy by instead minimizing the 
total thermodynamic potential 



HM,Nb) 



Fb{M,Nb)- 
+Ils{Co)M 



(11) 



over M and A''^ with Cq held fixed. This is analogous 
to minimizing the grand potential F — fiN, but allows 
for variation of the polymer layer's volume as well. 
By direct differentiation it can be verified that the ex- 
trema of $ correspond to equal chemical potential and 
equal osmotic pressure in the layer and in solution. 



B. Free Energy Landscape 

Setting /is and fis equal yields an explicit solution 
for 0. This gives, in equilibrium. 



(1-^)7 
7 + exp(xV') 



(12) 



with7 = Co/(l-Co). 

If one makes use of Eq. ( 12 1 to reduce <I>(Af, Nb) to 
one variable, one obtains 

$ = h'^/{2N) + M{l-ij)\n[l-ij] 

-Mln[l-Co] (13) 
-M(l - 0;) ln[l +7exp(-x^)]. 

Note that the first two terms are simply the free en- 
ergy obtained for a polymer brush with no nanoparti- 
cles, Eq. Q. Since = N/M, this form allows $ to be 
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regarded as a function of one variable, M (or, equiv- 
alently, h = Mja^ where a is the distance between 
grafting sites), with Co as a fixed input parameter. 

It is instructive to plot $ vs. h at fixed Cq as seen 
in Fig. [ijfor X = —6 and x = ^ 10 (in units of kTb'^). 
For weak binding strength (e.g. x — ~6 as seen in 
Fig. |la[ ), $ vs. h for fixed Co has a single minimum 
for all Cq. As a result, the height of the polymer layer 
changes smoothly upon the addition of more nanopar- 
ticles. On the other hand, if the binding strength |x| 
between the nanoparticles and the chains exceeds a 
certain threshold, the free energy can have two local 
minima as seen in Fig. |lb[ with the global minimum 
representing the equilibrium state. Within mean field 
theory, this results in a discontinuous transition in the 
polymer layer between a high h (extended) and low h 
(collapsed) state when sweeping through Co- 





Figure 1: The total thermodynamic potential $ vs. 
polymer layer height h ioT a = A, N = 100, and (a) 
X = -6 and lO'^ < Co < 10~^ or (b) x = -10 and 
10^"* < Co < 10~^. In (a) there is a single minimum 
for all Cq. In (b) there exists a range of Co for which 
the potential has a double-minimum structure. At 
approximately Co ~ 1.2 x lO^'^ (shown via the bold 
red line), the two minima are equal in energy (i.e. a 
phase transition occurs). 



C. Behavior of the polymer layer's height as a 
function of particle concentration. 

In this subsection we show the predictions of the 
mean field theory of subsection |II A| regarding the ef- 
fect of the nanoparticles on the extension of the poly- 
mer brush layer at equilibrium. At moderate nanopar- 
ticle concentrations, the nanoparticles typically cause 
the brush to attain a more compact state. This de- 
crease in h may occur via a smooth decrease (for weak 
binding) or a first-order phase transition (for strong 
binding). Fig. [2] illustrates the effect of nanoparticles 
on the height of the polymer layer for various binding 
strengths. 

For sufficiently weak binding, such as x = —0.75, 
the effect of nanoparticle addition is only noticeable 
at very high nanoparticle concentrations. The rise in 
layer height at high concentration occurs because the 
brush must swell in order to accommodate the large 
number of bound nanoparticles. 

For moderate binding strength, such as x = —4, 
there is a range of Cq (approximately 10~^ < Cq < 
10^^) in which the height of the polymer layer de- 
creases smoothly as the nanopartcle concentration in- 
creases. Physically, if the brush is in a more compact 
state, each nanoparticle will interact, on average, with 
more monomers. These attractive interactions help to 
overcome the excluded volume repulsion between the 
polymers and thus favor the compact state. 

For sufficiently strong binding, such as x = 
—8, —13.5, this decrease in height occurs via a sharp 



transition, as illustrated in Fig. lb 



Note that a decrease in the height of the poly- 
mer layer upon addition of binding particles is qual- 
itatively consistent with the experiment of Ref. [T], 
where addition of the binding proteins causes a brush 
of biopolymers to strongly compress. 




Figure 2: Layer height, normalized to its value in the 
absence of nanoparticles /iq, as a function of the 
nanoparticle concentration in the bulk solution for 
several values of binding strength |x| with grafting 
distance a = 4 and N — 100. 

It might be thought that the nanoparticles simply 
play the role of solvent molecules, and their presence is 
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analogous to degrading the solvent quality. However, 
in a standard mean field analysis of a polymer brush in 
a variable solvent, one would start with the standard 
brush free energy in Eq. Q and add a term 
where x is the interaction parameter, to account for 
solvent quality. Such a free energy does yield a com- 
pression of the brush for poor solvent quality (x < 0), 
but does not contain the phase transition behavior 
predicted here. 



III. SIMULATIONS 

In order to augment and verify the mean field the- 
ory results, we now turn to coarse-grained Brownian 
dynamics simulations (e.g. [HJ of the system. 
In this type of simulation, a particle performs over- 
damped biased diffusion under the influence of the 
deterministic force arising from potential interactions 
with other particles and the random force that repre- 
sents the thermal fluctuations due to random collisions 
with the solvent molecules. In this paper, we are only 
interested in equilibrium properties, so that the sys- 
tem, comprised of monomers of the polymer chains 
and nanoparticles, samples the Boltzmann distribu- 
tion of conflgurations at the prescribed temperature. 
Hence the precise value of the friction constant is not 
critical, and we can assume sufficiently high friction so 
that an overdamped Brownian (Langcvin) algorithm 
is appropriate. 



A. Interactions 

We modeled the polymer strands as chains of beads 
with nearest-neighbor beads connected by finitely ex- 
tensible, nonlinear, elastic (FENE) springs |43], which 
exert a force of the form 



~kr 



FENE 



1 



(14) 



where k is the spring constant, r is the distance be- 
tween beads, and imax is the maximum separation of 
beads, at which the FENE restoring force becomes 
infinitely strong. This force behaves approximately 
as a Hooke's Law spring for r <C tmax, but becomes 
very large when r approaches Imax- As a result, it 
may be used when spring-like forces are expected to 
produce the correct physics, but large separations are 
physically unrealistic. 

In addition, all beads interact through a pairwise 
truncated 6-12 Lennard- Jones (LJ) potential [45 that 
represents short-range repulsion (excluded volume). 
We cut off the LJ potential at its minimum so that 
it is purely repulsive. Mathematically, the potential 
has the form 



Utrn.n. — 



In this potential, e sets the strength of the repulsive 
force, and b sets its spatial range (and thus the bead 
size). Since this LJ potential sets the monomer size, 
b is taken to be the same as the unit of length in the 
mean field theory. 

Nanoparticles are modeled as single beads of the 
same size as polymer beads. They, too, interact re- 
pulsively with each other via Eq. ([T5|), but they feel 



no analog of Eq. ( 14 ) and therefore do not form chains. 

We consider the case in which nanoparticles and 
polymer beads interact with each other attractively. 
This is achieved via a LJ-type potential which, in addi- 
tion to the repulsive part, has an attractive part com- 
posed of another LJ potential with a different value of 
e. The two parts are pieced together continuously at 
their common minimum (r = &). That is. 




£6, 



r <b 
r > b 



(16) 



where et is the binding energy of the attractive 
nanoparticle-monomer interaction. The potential in 
Eq. ( 16 1 ensures that the strength of the attractive in- 



0, 



> b 



(15) 



teraction can be varied without simultaneously chang- 
ing the excluded volume interaction. 

Nanoparticles and monomers both interact with the 
grafting surface via another truncated, purely repul- 
sive LJ-type potential centered at ^ = and truncated 
at z = b. This force depends on the z coordinate of 
the particles only, and so models a flat wall. The po- 
sitions of the terminal beads of all chains are fixed 
at z = 0, so that each chain is rigidly grafted to the 
surface at one end. 

The results reported here are for chains of length 
A'^ = 100 beads that were grafted onto a square lat- 
tice with lattice parameter a = 4 or a = 3, giving a 
grafting density of cr = 1/a^. Periodic boundary con- 
ditions were used in the x and y directions. Unless 
otherwise noted, a system size of 16 chains (4 x 4) was 
used. 



B. Simulation Measurements and Procedures 

The simulation box contains the grafted polymer 
layer, above which there is a solution of nanoparticles 
free of polymers. The upper boundary ("ceiling") of 
the simulation box is placed in contact with a reservoir 
of nanoparticles. The solution of nanoparticles above 
the brush is held at fixed concentration (within fluc- 
tuations) by flxing the incoming flux of nanoparticles 
from the reservoir, while the outgoing flux is deter- 
mined by making the "ceiling" an absorbing bound- 
ary. This means that the entire simulation box comes 
to equilibrium with a reservoir of nanoparticles at a 
fixed chemical potential (i.e. a grand canonical en- 
semble is simulated). For fixed flux into the box, it 
was verified that the time-averaged concentration re- 
mained constant after equilibrium was reached, and 
this concentration was found to be independent of the 
properties of the brush (e.g. h, et). 
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The brush height was measured by creating his- 
tograms of the monomer density as a function of z 



such as those in subsection IV C and marking the top 
of the brush as the z value at which the monomer 
density became neghgible. Once h is determined from 
the monomer density profile, all nanoparticles with 
z < h are considered bound, and those with z > h are 
considered to be in solution. 



IV. SIMULATION RESULTS 

A. Behavior of the polymer layer height as a 
function of particle concentration 

Simulation results support the mean field picture, 
which was illustrated in Fig. [2] Using the same graft- 
ing density as in Fig. [2] as well as a second grafting 
density, a — 3, the simulations produced data shown 
in Figs. |3]and|4j Fig. [Sjoverlays the simulation results 
and the mean field theoretical results shown in Fig. [2] 
(using the same values of x)- One fitting parameter 
has been used for each €b to match the MFT and sim- 
ulation results, namely the x value that corresponds 
to each ej,. x should only depend on the interaction 
strength eb and the volume around each nanoparti- 
cle in which interactions occur, not on other details 
of the system. Hence, the same x fitting parameters 
were used in both Figs. [3] and |4] demonstrating that 
for a given e^, x does not depend on a. 
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Figure 3: h/ho vs. Cq for four different binding 
strengths and grafting distance a — -i. Grand 
canonical simulation results (circles) and MFT 
results (lines) plotted on the same axes. The x 
values for the analytical lines are the same as those 
of Fig. [2] and correspond to the values shown in 
the legend. 

Fig. [5] shows the nanoparticle concentration in the 
polymer layer </> vs. the nanoparticle concentration 
in solution Cq. This figure also shows qualitative 
agreement between MFT and simulation, including 
the presence of a phase transition for = 2 and 
X = —13.5. The quantitative discrepancies between 
MFT and simulation are largest when the number 
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Figure 4: h/ho vs. Cq for four different binding 
strengths and grafting distance a = 3. Grand 
canonical simulation results (circles) and MFT 
results (lines) plotted on the same axes. The x 
values for the analytical lines are the same as those 
of Fig. [2] and correspond to the eb values shown in 
the legend. 



of nanoparticles is small and the binding strength is 
high. In this limit, there will be strong correlations be- 
tween the locations of nanoparticles and monomers, 
and hence the mean-field picture of uniform density 
and random mixing begins to break down. 




Figure 5: (j> vs. Cq for four different binding 
strengths and grafting distance a = 4. Grand 
canonical simulation results (circles) and MFT 
results (lines) plotted on the same axes. The x 
values for the analytical lines are the same as those 
of Fig. [2] and correspond to the values shown in 
the legend. 

From the fitting of a % value to each Cf,, we found 
that Ixl increases approximately linearly with Cfc. This 
is a sensible result since the free energy of interac- 
tion, a?hx4>'4' according to Eq. ([5|, should in simula- 
tion be comparable to the number of interacting bead 
pairs times the strength of the interaction, which is 
on the order of — e^. The number of interacting pairs 
is roughly the number of nanoparticles Nb = aP'hcj) 
times the number of monomer neighbors that each 



7 



nanoparticle interacts with. The number of monomer 
neighbors would be the volume around the "tagged" 
nanoparticle in which interactions occur, Vint, times 
the density of monomers ijj. This means that the free 
energy of interaction in simulation should be roughly 
—a?hti,Vint4>i^- Equating the mean field and simula- 
tion energies suggests x ~ —Vint^b- Since the effective 
range of the LJ force is of the order of 26, one expects 
Vint to be of the order of 2^ in units oi b^. A linear fit 
for xi^b) yields x = — 8.45efc + 4. Of course, this ar- 
gument is overly simplistic, so only rough agreement 
can be expected. One could also compare xi^b) to the 
second virial coefficient for the interaction. Although 
the dependence of the second virial coefficient on 
is nonlinear, the range of ef, considered here is small 
enough that the simplified linear approximation will 
suffice. 



B. Phase Transition 



the other hand, the curve for e?, = 1.5 lacks such a 
vertical drop, leading us to conclude that eb = 1.5 is 
not large enough for a phase transition to occur. This 
is consistent with the picture of a first order phase 
transition existing for sufficiently strong binding, but 
not for weak binding, as suggested by the mean field 
considerations. 
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For sufficiently high interaction strength, mean field 
theory predicts a discontinuity in the monomer and 
nanoparticle densities at the first-order transition. 
This prediction is consistent with the simulation re- 
sults of Figs. [3]and|4] 

In order to probe the region very close to the phase 
transition, it proved convenient to perform simula- 
tions with a fixed number of nanoparticles. In ad- 
dition to speeding up simulations and providing eas- 
ier access to situations in which very similar values of 
Co may yield very different equilibrium states, these 
simulations may be relevant to experimental setups 
in which the system size is small and the number of 
available nanoparticles in limited. We refer to these 
as "canonical ensemble simulations." In this type of 
simulation, nanoparticles are still free to partition be- 
tween the polymer layer and solution above, but they 
cannot enter or leave the simulation box. In this case, 
instead of being fixed initially by the boundary con- 
ditions, the concentration of free nanoparticles above 
the brush Co was determined after equilibrium was 
established by tracking the mean number of nanopar- 
ticles which partitioned into the solution. 

Fixing the number of nanoparticles enabled us 
to observe the layer morphologies in the transition 
regime in which the polymer layer could not (ac- 
cording to mean field theory) attain either the ex- 
tended phase or the collapsed phase globally. Analo- 
gies to simpler systems would suggest that this situa- 
tion could result in a coexistence region characterized 
by phase separation. 

As one can see in Fig. |6j canonical simulations and 
grand canonical simulations agree at concentrations 
for which data in available in both ensembles. In ad- 
dition, for £{, = 2 the canonical simulations attain the 
same value of Cq for a range of different values of Nb 
and h. This is a signature of a vertical drop in h vs 
Co, leading us to conclude that there is a phase tran- 
sition for eb = 2. The drop occurs at Co = Cq, which 
is estimated to be between 8 x 10~^ and 1 x 10"**. On 



10-' 10-* 10-' 10° 

Co 

Figure 6: Grand canonical simulation results from 
Fig. [3] overlaid with simulation results obtained by a 
canonical ensemble simulation method in which Co is 

not constrained by the boundary conditions. Both 
sets of data are for a = 4. 



Although most simulations were performed on a 
system of 16 polymer chains (4x4), the canonical sim- 
ulations were also performed with 36 chains (6 x 6) in 
order to better visualize the system and test for finite 
size effects. The 36 chain simulations (not shown) also 
show a vertical drop in h vs. Co with the same Cq as 
the 16 chain simulations (to within the uncertainty 
range indicated above). 

Snapshots taken from equilibrated simulations can 
be used to show morphology changes in the phase 
transition region. Fig. [7] shows snapshots taken from 
four of the simulations considered in Fig. |6] For a 
concentration just below Cq (Fig. [7^), the polymer 
layer appears to be homogeneous and extended, while 
for a concentration just above Cq (Fig. [t};), the layer 
appears to be homogeneous and collapsed. The canon- 
ical ensemble data allows us to access an intermediate 
case. Fig. [7)d shows a brush with eb — 2 which has a 
number of bound nanoparticles such that it has an in- 
termediate height {h/ho = 0.81). This snapshot does 
not show the system to be in a single, homogenous 
phase, but rather shows a region of high density at 
low z with a region of lower density above it. In con- 
trast, Fig. [7[i shows a snapshot of another polymer 
layer with a lower binding strength {eb = 1.5) and a 
number of nanoparticles such that h/ho — 0.84. In 
this case, the polymer layer appears to remain in a 
homogeneous phase despite its intermediate height. 
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Figure 7: (a)-(c): Snapshots from a simulation with 
ef, = 2 and a = 4 at various values of Cq. (a): Just 
below the phase transition (Cq — 6.9 x 10~^) the 
brush is extended and appears homogeneous, (b): At 
an intermediate brush height {h/ho — 0.81, 
Co = 1.0 X 10~^) the brush appears to display 
spatial inhomogeneities. (c): Just above the phase 
transition (Co — 2.4 x 10^"*) the brush is collapsed 
and appears homogeneous, (d): A snapshot from a 

simulation with ef, = 1.5 and a = 4 at an 
intermediate brush height {h/ho = 0.84) lacks the 
obvious spatial inhomogeneities of (b). 



C. Density Profiles 

When the nanoparticle concentration is small 
enough that they constitute only a small perturbation 
on the standard polymer brush, the monomer density 
profile is approximately parabolic as seen in Fig. [8] 
The height of the brush is approximately the inter- 
cept of the fit parabolia as shown in Fig. [Sj Although 
our MFT takes the monomer density to be constant 
throughout the polymer layer, the parabolic density 
profile obtained in simulations matches the expected 
behavior of a plane-grafted brush [T^], including 
the presence of a depletion region near the wall and a 
"foot" at high z in simulations [101 [TT] . 

Under certain conditions, the brush collapses, and 
space is almost completely filled with particles. As a 
result, the monomer density profile becomes approxi- 
mately a step function as seen in Fig. |9] Of course in 
simulations the step is not infinitely sharp, but instead 



includes a transition region of intermediate monomer 
concentration. We include this intermediate region as 
part of the brush. 




10 20 30 40 50 60 70 



Figure 8: A time averaged histogram of the 
monomers' height above the grafting surface for 
parameter values for which few nanoparticles are 
bound to the brush. The monomer density profile is 
nearly parabolic as expected for a brush with no 
nanoparticles. The red line shows a parabola to 
guide the eye. This figure was generated for a = 4, 
Ef, = 2, Co = 4.2 x 10-5. 



140 




5 10 15 20 25 

z 

Figure 9: A time averaged histogram of the 
monomers' height above the grafting surface for 
parameter values for which many nanoparticles are 

bound to the brush. The red line shows a step 
function which drops to zero at the measured value 
of h. The monomer density profile is close to a step 
function because monomers and nanoparticles are 
tightly packed as in a solid. This figure was 
generated for a = 4, ef, = 2, Co = 2.3 x 10^^. 

Typical nanoparticle density profiles corresponding 
to Figs. |8] and [9] can be seen in Figs. [TO] and [TT] In 
both cases, the nanoparticle density profile for z < h 
qualitatively tracks the monomer density profile be- 
cause nanoparticles in the polymer layer are expected 
to be bound to one or more monomers. 
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Figure 10: A time averaged histogram of the 
nanoparticles' height above the grafting surface for 
the same parameters as Fig. l8| 
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Figure 11: A time averaged histogram of the 
nanoparticles' height above the grafting surface for 
the same parameters as Fig. l9| 



V. SUMMARY AND DISCUSSION 



The results of mean field theory and Brownian dy- 
namics simulations not only show that the polymer 
layer is expected to attain a more compact state, they 
also reveal the possibility of a discontinuous collapse of 
a polymer brush in the presence of attracting nanopar- 
ticles as a function of nanoparticle concentration in 
the abutting solution. For sufficiently weak binding 
the collapse of the brush is smooth, while for suffi- 
ciently strong binding, the collapse occurs via a first- 
order phase transition from an extended state to a 
collapsed state. 

These observations may be useful in the field of 
polymer science. Control of polymer morphology in 
the presence of additives not only plays a role in the 
development of tougher polymeric materials (via the 
use of antiplasticizers, for example), but polymeric 
materials (sometimes decorated with nanoparticles of 
various kinds, e.g. [51 1161 are also becoming in- 
creasingly useful for applications such as organic elec- 
tronics. 

In addition, the insights gained in this study inform 
future analysis of the behavior of grafted biopolymers, 
such as those in the Nuclear Pore Complex (NPC) 
in which binding proteins play the role of nanoparti- 
cles [ini EO] ■ Of course, in physically realistic situa- 
tions, the theory and simulations presented here may 
be overly simplistic. In the relevant biophysical sys- 
tem of the nuclear pore complex, which is involved 
in the transport of proteins and nucleic acids into and 
out of the cell nucleus, the "nanoparticles" (which cor- 
respond to binding proteins) are significantly larger 
than those of our simulations [15], and they possess 
discrete binding sites rather than the isotropic bind- 
ing potential used here. Complications like these will 
be addressed in future work. In spite of its simplicity, 
the model presented here appears to share qualitative 
similarities with the NPC analogs, such as the collapse 
of nucleoporins upon binding to karyopherin transport 
proteins which was observed by Lim et al. [T]. 



In this paper, we considered the case of grafted 
polymers in equilibrium with nanoparticles that are 
attracted to the polymers. The mean field theory of 
and the Brownian dynamics simulations of 
are in qualitative agreement that, at certain 
interaction strengths and nanoparticle concentrations, 
the grafted polymer layer is expected to attain a more 
compact state upon the addition of binding nanopar- 
ticles. Hence, our results are in qualitative agreement 
with the biologically-motived experiments of Ref . [T] . 
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